Spin-orbit interactions mediated negative differential resistance in a 
quasi-two-dimensional electron gas with finite thickness 
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Effects of the spin-orbit interactions on the energy spectrum, Fermi surface and spin dynamics 
are studied in structural- and bulk-inversion asymmetric quasi-two-dimensional structures with a 
finite thickness in the presence of a parabolic transverse confining potential. One-particle quantum 
mechanical problem in the presence of an in-plane magnetic field is solved numerically exact. Inter- 
play of the spin-orbit interactions, orbital- and Zeeman-effects of the in-plane magnetic field yields 
a multi- valley subband structure, typical for realization of the Gunn effect. A possible Gunn-effect- 
mediated spin accumulation is discussed. 
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One of the major goals of modern electronics is the 
search for new spin-involved functionalities 1] in solid 
state nanoscale devices. Electric field controlled Rashba 
spin- orbit (SO) coupling which is originated from 
large potential gradient on the semiconductor/insulator 
interface of quantum wells and MOSFETs in the presence 
of macroscopic structural inversion asymmetry (SIA), is 
a promising tool [Sj in a realization of spin transport 
devices. On the other hand, special semiconduction ma- 
terials with a bulk inversion asymmetry (BIA) in their 
crystalline structure produce so-called Dresselhaus SO 
interaction Q , which interplays with Rashba interaction 
yielding unusual physical effects 0, @] . 

Different experimental techniques have been developed 
recently to control a coupling of spin to the electric field 
0-01 ■ An efficient g-tensor modulation resonance, ob- 
served in a parabolic Al^Gai^xAs quantum well 0] with 
varying Al content x = x{z) across the well, provided 
an opportunity to manipulate electron spins by means of 
various electron spin resonance type techniques. In-plane 
magnetic field in all of these experiments seems to be 
rather favorable to get a pronounced spin resonance. On 
the other hand, gate voltage control of spin dynamics in 
practical devices was shown |10|| to be reasonable only for 
a finite thickness of the electron gas. SO interactions in a 
2D electronic system produce an effective in-plane field, 
which results in an drift-driven in-plane spin polarization 
External in- plane magnetic field appears to be not 
always algebraically summed with SO induced effective 
field and results in a surprizing out-of-plane spin polar- 
ization [l^ , which was observed in strained n — InGaAs 
film 's'l . On the other hand Hanle precession of optically 
oriented 2D electrons in GaAs can be well described 
by algebraic addition of these fields. All these facts show 
nontrivial effects of in-plane magnetic field on spin dy- 
namics in quasi-2D systems. 

In this paper we show that an interplay of Rashba and 
Dresselhaus SO interactions in a quasi-two-dimensional 
(quasi-2D) electron gas with finite thickness results in a 



multi-valley ('N'-shaped) energy dispersion in the pres- 
ence of an in-plane magnetic field. The obtained energy- 
momentum relation looks like those existing in GaAs 
and InP-type materials of the transferred-electron de- 
vices, where negative differential resistance and Gunn mi- 
crowave current oscillations take place. Finite thickness 
of the electron gas gives rise to transverse-quantized lev- 
els with Rashba-Dresselhaus spin-split subbands in each 
level. The in-plane magnetic field removes degeneration 
in the anticrossing points of different subbands, open- 
ing a magnetic field-controlled gap at these points. Bar- 
rier height between the valleys in the 'N'-shaped energy- 
momentum relation and the valley curvatures seem to be 
controlled by Rashba- and Dresselhaus SO coupling con- 
stants, the in-plane magnetic field strength and Lande 
factor. The dependence of the energy spectrum on 
the magnetic field shows strongly non-linear and non- 
monotonic behavior even in weak field strengths. 

The single particle Hamiltonian of the system in the 
effective mass approximation is written as 



2^2 



H 



m ujnZ 



2m* 



eEgZ - 



1 



g^J,B(TB, (1) 



where P = p — -A is an electron momentum in the 

^ c 

presence of a vector-potential A, m* and e are the elec- 
tronic effective mass and charge, respectively; Eg is a 
strength of the gate electric field. A parabolic poten- 
tial with a frequency wq, characterizing the electron gas 
thickness, does not produce SO interactions, and is cho- 
sen in Eq.([T]) to confine the electron gas in z-direction. 
Since Rashba SO interaction in the conduction band of 
a semiconductor is determined by the electric field in the 
valence band rather than by that in the conduction band 
14| , the parabolic confinement approximation neglects a 
small interface contribution to Rashba SO coupling con- 
stant. The last term in Eq.([T]) is Zeeman splitting energy 
in the external magnetic field B with hujz ~ gpLBB/2, 
where fiB — the Bohr magneton of a free elec- 

tron with mass toq, g is the effective Lande factor, and 



c = WxjO'yjaz} are the Pauli spin matrices. Hso in 
Eq. ([Ij contains Rashba term as well as Dresselhaus term 
with characteristic parameters a and /?, correspondingly 
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The electronic wave function '^{x,y,z) can be ex- 

ik^x+ikyy (i>^(z)\ r 



pressed in the form 'i'(a::,2/, z) = 
a magnetic field aligned along x-axes B = {i?, 0,0} 
under the gauge A — {0,— i?2;,0}. In the absence 
of the SO interactions and Zeeman term, Hamil- 
tonian becomes diagonal and the equations for 
-01 and ipi are reduced to the oscillator equation 
i-LQipn\z) — En'ipri^ with real wave function ipn\z) = 
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and the energy spectrum E,, 
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Hn{z) is the Hermite polynomial. General solutions are 



chosen as linear combinations of ^/'i''^(z) 



E 



2^0 



as 



(3) 



where cr =t,i- By putting Eq.® into Schrodinger equa- 
tion one gets the matrix equations Ncra*^ — for the vec- 
tors sl" = {aQ,aJ,a2, . . . }. with a =tii are square 
penthadiagonal matrices of infinite order with non-zero 
entries ^ only if \i - j\ < 2, and Nt = (N^)*- 
The energy spectrum has to be found from the secular 
equation, by equating the determinant of the matrix N 
to zero. The infinite penthadiagonal matrix is truncated 
down to the first n rows and n columns, and the roots of 
its determinant are found by numeric methods. All pa- 
rameters in our numeric calculations, shown below with 
tilde, are done dimensionless in the unit of the charac- 
teristic frequency wo of the confining potential or in the 
length scale Iq = (S/m*wo)^^^, which is a measure of the 
electron gas thickness. 

Energy spectrum. The solutions of the n x n determi- 
nant converge very rapidly as n increases. The energy 
dispersion for the first three transverse quantized levels 
(n=3) is depicted in FiglT] for different values of the SO 
coupling constants and the magnetic field. In the absence 
of the in-plane magnetic field and one of the SO coupling 
constant, the energy spectrum En{ky) is described by the 
well-known two symmetric parabolas in each transverse- 
quantized level (dashed (green) curves in Fig[T^). An 
external gate electric field fg = eEgl^/ {ojqK) coherently 
shifts the energy subbands. Non-zero kx component of 
the momentum splits spin-up and spin-down spectra at 
fcy = 0. The energy spectrum becomes asymmetric along 
momentum- and energy-axes if q = a/(fia;o'o) 7^ 0, 
/? = f3/{hujolo) 7^ and kx 7^ 0. Energy spectra cor- 
responding to the same transverse-quantized level with 
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FIG. 1: (Color online) Dependence of the energy spectrum 
on Rashba- and Dresselhaus SO coupling constants as well 
as on the orbital- and spin-effects of the in-plane magnetic 
field. E vs. ky is shown (a) in the absence of the magnetic 
field at kx = 0, fg = 0.3, a = 0.5, P — hy dashed (green) 
curves; kx = 0.4, fg — 0, a = 13 = 0.6 by dot-dashed (blue) 
curves and kx = 0.4, fg — 0, a = 0.8, /3 = 0.4 by solid 
(red) curves; (b) under the in-plane magnetic field ujb ~ 0.7, 
g — 0.2 and kx = 0.2, /3 = 0.4 but for different values of a 
and fg: a = 0.4, /g = by dashed (green) curves, a — 1.2, 
fg — 0.3 by dot-dashed (blue), and a = 0.9, fg — hy solid 
(red) curves. 

opposite spin orientations split at the intersection point 
for a ^ whereas they do not split at the intersection 
point for a = 13. Nevertheless, anticrossing of the energy 
spectra corresponding to different levels persists irrespec- 
tive of the values of a and /?. FiglT}) describes the depen- 
dence of En on ky in the presence of the in-plane magnetic 
field ujb = wb/wq, Zeeman splitting g = gm* / (Amo) and 
kx for different values of the SO coupling constants. 
The magnetic field seems to remove the degeneracy at 
the anticrossing points, and to open a gap, which in- 
creases with magnetic field. Although the Zeeman split- 
ting strongly modifies the energy levels in the presence 
of only Rashba SO coupling, shifting the parabolas bot- 
toms along ky and splitting the energies of spin-up and 
spin-down electrons, it does not change the levels sym- 
metry in the presence of only Dresselhaus SO coupling 
(dot-dashed (blue) curves in FiglT})). The same picture 
is obtained for the energy dispersion along kx under the 
similar conditions given in FigH}) but with the replace- 
ments kx -H- ky and a O /3. 

The dependence of the energy spectrum on the in-plane 
magnetic field under different momentum components 
and SO coupling constants is given in Figl2}i for zero 
Lande factor and in Figl2}) for a non-zero Lande factor, 
g — 0.2. The energy spectrum in both cases displays 
considerably nonlinear and non-monotonic behavior even 
under weak magnetic fields, when < ujb < 1. The mag- 
netic field splits the spectrum at the anticrossing points. 
The energy dispersion in the negative magnetic field is 
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FIG. 2: (Color online) Dependence of the energy spectrum 
on the in-plane magnetic field (a) in the absence g = and 
(b) in the presence g = 0.2 of Zeeman splitting at fg — 0. 
The parameters kx — 0.4, ky = 0.5, a — 0.6, and /3 = 0.4 
correspond to dotted (green) curve. Dashed (blue) and solid 
(red) curves depict the cases kx = 0.2, ky — 1.4, a = 0.8, 
P = 0.4 and k^ = 0.4, ky = 1.8, a = 0.6, (3 = 0.4 respectively. 

exactly the same as the dispersion in the positive mag- 
netic field but with reversed in sign the Lande factor. 

Fermi surface. Influence of different SO interactions 
and in-plane magnetic field on the Fermi surface is 
shown in Fig|3]for electrons, resided the first transverse- 
quantized level. The dimensionless energy is fixed in 
the middle of the two adjacent transverse levels Ep = 
Ep/{hujQ) = 1.0. In the presence of only one SO in- 
teraction the Fermi surface is symmetrically splitted at 
ujb = into two surfaces of spin- up and spin-down elec- 
trons even under non-zero gate voltage. Fermi surfaces 
become strongly anisotropic for a ^ and P ^ due 
to shifting of two circles away from each other in both 
kx and ky directions. If a ^ ^, the circles correspond- 
ing to spin-up and spin-down electrons split at the in- 
tersection points. Nevertheless they do not split at the 
intersection points if Rashba and Dresselhaus constants 
are equal each other, a = (3. Gate voltage increases the 
anisotropy, but it does not split the Fermi surface at the 
intersection points for a = (3. Instead, in-plane magnetic 
field removes a degeneracy at the intersection points. 

Spin-Gunn effect. Multi- valley energy dispersion, dis- 
cussed above, allow us to suggest an existence of neg- 
ative differential resistance and spin-Gunn effect [3] in 
a quasi-2D electron gas in the presence of SO interac- 
tions. Assume that the Fermi level crosses the second 
transverse-quantized level in the energy-momentum dis- 
persion, shown in FiglTJj with solid (red) curve under the 
conditions h = 0.2, & = 0.9, /3 = 0.4, wj? = B/Bq = 0.7, 
and g = 0.2. In order to estimate a, /3, Bq — m*cujo/e 
and the electron gas thickness Iq — y^h/m*ujQ, we use 
a ~ 1 X 10~^eV ■ cm, and /? ~ 0.4 x lO'^eV ■ cm, as 
well as ujQh ~ 10 meV, [lB|, which yield Iq = 10.6 nm 



FIG. 3: (Color online) Fermi surface of electrons in the first 
transverse-quantized level, (a) Dot-dashed (green) curves cor- 
respond to Q = 0.5, P — 0, fg — 0.3, dashed (blue) curves to 
a = 0.5, /3 = 0.4, fg~0 and solid (red) curves to a = (3 — 0.5, 
/g = in the absence of magnetic field lob = g — 0; (b) Dot- 
dashed (green) curves correspond to Cjb = § = 0, fg = 0.3 for 
a = P = 0.5, dashed (blue) curves to ljb = 0.7, g = 0, fg — 
for a — P — 0.5, and solid (red) curves to Cjb ~ 0.7, g = 0.3, 
/g = for a = 0.5, jS = 0.4. 

in agreement with quantum well width in the experi- 
ments 0, 13, Bo = 0.9 T and B = 0.6 T (correspond- 
ing to UB = 0.7), a = 0.944 and /3 = 0.38 for, e.g. 
conduction band electron in GaAs with m* = O.OGSmg. 
Indeed, Rashba SO constant varying in the interval of 
a ~ (4.47 — 6.30) x 10~® eV ■ cm has been reported 17| 
for InAlAs/InGaAs heterostructures. Furthemore a gi- 
ant SO splitting was recently observed for quantum well 
states of a Bi monolayer on Ag(lll) fl^, on Si(lll) |l9j, 
and on Cu(OOl) '^2(1] with Rashba parameters correspond- 
ingly a = 3.05 x 10'^ eV ■ cm, 1.37 x 10"^ eV ■ cm and 
(1.5-2.5) X 10-^ eV -cm. 

In equilibrium most electrons reside near the bottom 
Eki of the lower valley in Fig|3^, which corresponds to 
an anticrossing point with small effective mass m^ and 
higher mobility /ui due to the sharp curvature of the 
valley. The electric field, applied along ky axes, bends 
the Fermi level and accelerates the electrons with a def- 
inite (up- or down-) spin polarization to the top of the 
valley separation Em when the field reaches a thresh- 
old value. The upper valley is more fiat where the ef- 
fective mass 7712 is rather heavier and the mobility 112 
is smaller than /ii. The ratio 7712 /?7i J is estimated to 
be ~ 80 by fitting the valleys' bottoms to parabolas. 
Our calculations show that the effective masse in a val- 
ley, obtained by anticrossing of a parabola of 71th sub- 
band with a parabola of (71 -|- l)th subband, can be con- 
trolled. Indeed, by changing the signs or the amounts 
of SO coupling constants, magnetic field, gate voltage 
and g-factor the bottoms of the parabolas are control- 
lably shifted in different directions along the energy- and 
momentum-axces, giving different curvatures and effec- 
tive masses in the valleys. The total electron density 77 is 
written as a sum of the densities ni and 712 in the lower 
and upper valleys, correspondingly, n — ni + n2. The 
population ratio between the upper and lower valleys is 
assumed to be given by Maxwellian energy distribution 
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FIG. 4: (Color online) (a) Two valley spectrum correspond- 
ing to the lower branch of the first and the second transverse- 
quantized levels with a — 0.9, /3 — 0.4, ub = 0.7, g — 0.2, 
fg~0 and kx = 0.2 shown by solid (red) curves, and a = 1.2, 
P = 0.4, CbB = 0.8, g = 0.3, fg = 0.3 and L = 0.4 shown by 
dot-dashed (blue) curve. The dependence of (b) the average 
drift velocity v and (d) the spin accumulation P on the electric 
field £ for /ii = 10 m'^/V-s and T = 4.2K, = /.12/M1 = 0.02, 
R = 50. (c) P vs. £ for m = 20 m'^/V • s and T = 50K, 
fi = ^12/ ^ii = 0.1, R = 10. Four different values of the activa- 
tion energy /S.E = 4 meV, 6 meV, 8 meV and 10 meV corre- 
spond to dotted (green), dashed (violet), dot-dashed (blue), 
and solid (red) curves in each figures. 

as n^jrii — Reyj^{~ jkT^^ where Te is the electron 
temperature, i? is the ratio of the density of states in 
the upper A'2 and lower A^i valleys with i? — m\lm\ for 
quasi-2D electron gas, and Ai? — Ek2 — Eki is the energy 
separation between the two valley minima, which can be 
estimated from e.g. FigH^ to be AE = 0.45 fkuQ. Then, 
electron accumulation P with definite spin polarization 
in the upper valley is written as 



„ no 
P = — = 



n2 



Rejip{-AE/kTe) 
ni+n2 ^ 1 + Rexp{-AE/kTe 



(4) 



The electron temperature is defined from the condition 
eSvT^ = 3k{Tf, — T)/2, meaning that electrons gain an 
energy eSvTe from the electric field £, aligned in y-axis, 
as they move up the potential barrier with the average 
drift velocity v = (/iini -\- ^i2n2)£ / {ni + n2)^ and transfer 
the excess energy to the atomic lattice with temperature 
T in the energy relaxation time Tg, 2l|. The steady- 
state current J is expressed through v as J = e(/xirti -|- 
fi2n2)£ = env, which is calculated from the following 
self-consistent equations 



n = T- 



3^' 



l + i^i2/t^i)ReM-^E/kT,) 
1 + Rexp{~AE/kTe) 

(5) 

Eqs.(|4]) and ([5]) determine completely v — £ 01 I — V 
dependence, which displays a shape with a region of 



negative differential resistance. Numerical solution of 
Eqs.(|4]) and (O yields a dependence of the spin-polarized 
electron accumulation P on £, which is presented in 
FigHt, and d for temperatures T — 50K and A.2K, 
correspondingly, v vs. £ and P vs. £ dependences 
are depicted in each figures for different valley separa- 
tions AE = 4 meV, 6 meV, 8 meV and 10 meV, which 
are shown by dotted (green), dashed (violet), dot-dashed 
(blue), and solid (red) curves, respectively. High tem- 
perature mobility of electrons in, e.g. GaAs varies in the 
interval of (10-^ 20) m^/V ■ s [211 . The energy relaxation 
time Te is chosen to be = 10~^^ s. The higher valley is 
basically resided at T = 50K by electrons, according to 
FigHt, due to thermal fluctuations even in the absence of 
the electric field. The spin-polarized electron accumula- 
tion slightly increases from P ~ 0.6 up to the saturation 
value P ^ 0.90 as £ increases and reaches ^ 1.5 kV/cm. 
The spin-dependent electron accumulation at T = 4.2if, 
given in FigHJl for fii — 10 rr? jV • s, 112I [i-i = 0.02 and 
m\lm\ — 50, sharply increases with electric field from 
P 0.0 up to P - 0.98 at f 1 kV/cm. Fermi energy 
moves up with increasing the electron concentration and 
it crosses a subband with opposite spin polarization. 

In conclusion, we show that interplay of SO interac- 
tions and in-plane magnetic field strongly changes an 
electron spectrum of quasi-2D electron gas with finite 
thickness yielding multi-valley energy dispersion in each 
transverse-quantized subband. As a result, an instable 
regime of negative differential conductance is realized for 
spin-polarized electrons. 
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